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SYNOPSIS 


To  model  the  nuclear  induced  lightning,  the  basic  processes  responsible 
for  this  phenomenon  must  be  understood.  The  air  chemistry  plays  an  essential 
role  in  such  a  modelling.  However,  a  self  consistent  approach  to  the  modelling 
requires  the  incorporation  of  the  best  available  reaction  rates  and  cross  sections 
for  the  basic  processes.  A  review  of  these  processes  and  the  relevant  rates  is 
undertaken  for  the  sake  of  a  realistic  modelling  of  NIL. 


A  REVIEW  OF  THE  AIR  CHEMISTRY  AND  RELEVANT  PARAMETERS  FOR  THE 
MODELLING  OF  THE  NUCLEAR  INDUCED  LIGHTNING 

1.  INTRODUCTION 

Nuclear  induced  lightning  (NIL)  has  been  observed*  from  detonations  at 

the  sea  level  (sunface  bursts).  These  discharges  occur  some  distance  away 

(900  -  1400  meters)  from  the  burst  point  and  they  grow,  from  the  top  of 

structures  which  are  above  ground  by  some  meters,  upwards  to  the  cloud 

level.  The  growth  velocity  of  the  discharge  is  found*  to  be 

; 

~  1  -  2  X  10  cm/sec  and  that  the  discharge  luminosity  lasts  for  times 
of  -  0.1  sec. 

The  mechanism  for  the  discharge  propagation,  i.e.,  the  streamer 
initiation  and  propagation,  requires  the  presence  of  charges  and  electric 
fields.  These  physical  quantities  arise  from  the  interaction  of  y  -rays, 
produced  by  the  nuclear  detonation,  with  the  air  molecules.  The  interaction 
produces  Compton  electrons  which  move  outward,  ionizing  air  species  in  its 
path,  resulting  in  a  plasma  current  moving  in  the  opposite  direction.  These 
processes  are  understood  to  generate,  the  charges  and  electric  fields  which 
are  responsible  for  the  discharge  initiation  near  sharp  objects.  A  simple  and 
semi -quantitative  analysis  for  the  electric  field  has  been  performed  in  Ref. 
1,  where  the  the  radial  electric  field  is  estimated  to  be  300  Vcm“*.  This 
value  for  the  field  Is  much  below  the  air  breakdown  threshold  field^  of  30 
Vcm~*  Torr“*. 

The  NIL  phenomenon  is  complex  and  requires  a  multidisciplinary  approach 
for  Its  understanding  and  modelling.  The  modelling  should  comprise  the 
following  areas  of  physics  and  chemistry:  interaction  of  radiation  with  air 
species,  Ionization  of  air  and  generation  of  air  conductivity,  electromagnetic 
equations  to  solve  for  the  electric  field  in  conjunction  with  the 
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conductivity,  the  air  chemistry,  air  heating,  hydrodynamic  and  flow  equations, 
breakdown,  streamer  propagation,  and  discharge  radiation. 

A  comprehensive  modelling  of  NIL  in  time  and  space,  which  encompasses  all 
physics  aspects  discussed  above,  is  an  approach  which  could  lead  to  a  complete 
understanding  of  the  phenomenon,  even  though  such  an  approach  may  be  an 
ambitious  one.  A  step  in  this  direction  has  been  undertaken  by  Longmire,  et 
al^.  The  success  of  the  modelling  of  NIL  or  any  other  phenomenon,  however, 
depneds  on  a  large  measure  on  the  assumptions  made,  the  availability  and 

accuracy  of  the  relevant  data  used  for  the  basic  processes  involved.  For 
nuclear  induced  lightning  the  air  chemistry  plays  an  important  role.  This 
role  is  apparent  in  two  important  areas  of  the  lightning  discharge.  The 
first,  clearly  is  exhibited  in  the  y  -ray  generated  conductivity  in  the 
disturbed  region  of  the  atmosphere,  and  the  other  is  the  air  chemistry  at  the 
tip  of  the  sharp  object  from  which  the  streamer  propagates. 

This  report  deals  with  the  air  chemistry  relevant  to  NIL,  where  it 
reviews  the  rate  coefficients  for  the  apporpriate  processes  and  provides 
current  and  best  values  for  modelling  purposes.  Various  atmospheric  models 
exists  for  EMP  calculations  and  a  study^  has  been  made  to  determine  the 

uncertainties  In  these  calculations  due  to  various  chemistry  models  and 
coefficients.  Scheibe®  have  considered  the  generation  of  smog!  and  its  role 

in  the  electron  attachment  In  an  EMP  environment  at  the  sea  level.  Grover  and 

Gilmore®  have  reviewed  the  mobility  and  attachment  data  for  the  EMP 
environment. 

An  important  element  in  the  NIL  modelling  is  the  air  conductivity  whose 
calculation  depends  heavily  on  the  air  chemistry.  The  conductivity  is  a 
function  of  the  electron  density,  and  ion  density  when  the  conductivity  is  due 
to  Ions,  and  their  collision  frequencies  with  the  neutral  species.  The 
electron  and  ion  density  calculations  require  a  set  of  rate  equations  which 


determine  their  production  and  loss  mechanisms.  These  processes  are  numerous 
and  involve  a  large  number  of  reaction  and  their  coefficients  which  are 

temperature  dependent.  The  collision  frequency,  on  the  other  hand,  can  be 
obtained  if  the  appropriate  velocity  distributions  are  known  along  with  the 
relevant  cross  sections.  This  approach  requires  one  to  solve  the  Boltzmann 
equation.  Other  approaches  are  also  available,  e.g.  assuming  a  Maxwellian 
velocity  distribution  for  the  electrons,  or  utilizing  the  appropriate  particle 
mobilities  obtained  experimentally  as  a  function  of  E/p  and  the  average 
electron  energies,  also  obtained  experimentally  as  a  function  of  E/p.  The  air 
chemistry  near  the  tip  requires  a  detailed  treatment  of  air  heating  processes 
and  the  energy  flow  into  various  modes  of  freedom,  e.g.  kinetic,  vibrational, 
etc.  as  well  as  the  ionization  frequency  as  a  function  of  E/p. 

This  report  deals  with  the  prospect  of  a  comprehensive  NIL  modelling  and 
tries  to  provide  the  basic  data  for  the  most  important  processes.  These 
processes  are:  attachment,  ion-ion  recombination,  ionization  frequency, 
association  reactions,  dissociative  recombination,  vibrational  excitation  and 
relaxation,  the  electron  temperature,  thermal  processes  and  others. 

2.  ATTACHMENT 

The  electron  attachment  to  Og  in  air  proceeds  mainly  by  the  following 
processes 


e  +  02  +  0"  +  0 


e  +  O2  +  M  +  O2  +  M 


The  first  process,  the  dissociative  attachment,  has  a  threshold  of  -  4.3  eV 
and  hence  it  is  important  when  the  electron  energy  is  high.  The  cross 
section^  for  the  dissociative  attachment  has  been  averaged®  with  the  electron 
velocity  over  a  Maxwellian  electron  velocity  distribution  and  the 
corresponding  rate  coefficient  is  given  in  Table  1  as  a  function  of  the 
electron  temperature. 


Table  I  The  Dissociative  Attachment  Rate  Coefficient  for  0 


The  data  of  Table  I,  however,  can  be  expressed  (within  20%)  and  in  the 
temperature  range  of  (0-5)  eV  by  the  following  expression 

a2  =  1.0  x  10'11  [2.28  +  11.6  Tg  -  1.75  Tg2]  Te'0,5  Exp(-4.2/Te)  (3) 

Where  Tg  is  in  units  of  eV. 

The  three-body  attachment  process  as  given  by  Eq.  2  indicates  that  M  is 
the  third  body  which  may  be  O2,  N2,  H2O,  CO2  and  any  other  molecular  compound 
generated  in  the  disturbed  air. 

Considerable  measurements  on  the  three-body  attachment  rate  coefficient 
exists,  indicating  the  predominance  of  this  attachment  for  electrons  with 
energy9  of  <  1.0  eV  ..  The  Measurements,  prior  to  1971,  have  been  reviewed  by 
Phelps9’*9  and  recommended  rate  coefficients  for  attachment  to  O2  are  given 
for  various  third  bodies  appearing  in  Eq.  2.  These  and  more  recent 

measurements  which  are  reviewed  here,  have  been  performed  for  thermal 

electrons  and  over  a  limited  range  of  E/p.  And  are  mainly  obtained  from  swarm 

and  electron  beam  experiments. 

2.1  THREE  BODY  ATTACHMENT  (M=02) 

For  the  three  body  attachment  in  pure  O2  and  for  thermal  electrons  at 
300°K,  the  measured  values  of  various  workers  are  in  good  agreement.  These 
values,  in  units  of  cm6/sec,  are:  (2.8  +  0.5)  x  10"39  [Ref.  11],  2.0  x  10"39 
[Ref. 12],  (2.4  ±  0.1)  x  10-30  [Ref.  13  and  14],  (2.1  ±  0.2)  x  10'30  [Ref.  15], 
3.8  x  10’30  [Ref.  16],  2.1  x  10_3°  [Ref.  17],  -  2.5  x  10_3°  [Ref.  18].  It 

seems  that  a  value  of  2.4  x  10"39  cm®/sec  is  a  resonable  value  for  thermal 
electrons  at  gas  temperature  of  300°K.  The  temperature  dependence  of  the 


three-body  rate  coefficient,  where  O2  is  the  third  body,  have  been 

obtained  in  swarm  experiments  from  the  attachment  coefficients  per  unit 
length,  n  ,  which  are  measured^’^’^  as  a  function  of  jj-  (or  E/p)  .  In 
general,  83(02)  is  related  to  n  according  to  (4) 


n  =  a3 (02)  02 


where  Vd  is  the  electron  drift  velocity  which  is  a  function  of  E/N  (or  E/p) 

on 

and  has  been  measured  in  a  large  number  of  gaseous  media  .  The  attachment 
coefficient  per  unit  length  signifies  the  total  attachment  in  O2,  i.e.,  the 
two-body  and  the  three-body  processes  delineated  experimentally  where  the  two 
body  is  independent  of  pressure  and  important  at  high  E/p,  in  contrast  the 
three-body  is  pressure  dependent  and  predominates  in  the  regions  of  low  E/p. 
These  characteristics  can  be  seen  in  Figure  1  which  is  based  on  data  from 
Reference  19. 

The  three-body  attachment  rate  coefficient,  83(02)  can  be  obtained  from 
data  such  as  the  one  given  in  Figure  1,  Equation  3  and  the  drift  velocity  as  a 
function  of  E/p.  In  addition  one  needs  the  data  for  the  mean  electron  energy 
or  the  characteristic  energy  (~Te)  which  are  also  measured  as  a  function  of 
E/p.  From  these  considerations  one  obtains  83(02)  as  a  function  of  Te,  where 
most  of  the  data  are  basically  for  electron  temperatures  of  leV  and  lower  with 
the  gas  temperature  at  300°K.  However,  several  measurements  have  been  made 
where  the  gas  temperature,  Tg,  is  varied.  Such  data  are  as  follows:  Chanin, 
et  al  ^  for  Tg  between  77°K  and  370°K,  Pack  and  Phelps^  for  Tg  between  300°K 
and  575°K  and  Truby^  for  Tn  between  113°K  and  300°K.  Phelps^  recommends  the 


following  expression  for  the  three-body  attachment  for  thermal  electrons  when 
the  gas  temperature  is  varied  from  300°K  to  600°k 

a3  (02)  =  (1.4  ±  0.2)  x  10‘29  Exp  ^ 

On  the  other  hand,  the  measured  rate  coefficient  of  Truby10  fits  the  following 
expression . 

a3  (02 )  =  4  x  10‘30  Exp  (-  193/T  )  (6) 

However,  eventhough  the  measurements  of  Truby10  are  for  temperatures  below 
300,  his  expression  (Eq.  6)  can  be  used  for  higher  temperatures  (up  to  600°K) 
where  it  differs  in  value  by  10%  from  that  of  Phelps10  (Eq.  5). 

The  attachment  rate  coefficients  as  a  function  of  the  electron 
temperature  are  shown  in  Figure  2  based  on  various  swarm  data.  These  data 
show  a  smooth  variation  of  the  coefficient  with  the  electron  temperature. 
However,  the  data  from  the  electron  beam  experiment10,  also  shown  in  Figure  2, 
show  definite  structures  as  a  function  of  the  electron  temperature.  The  peaks 
of  the  structures  occur  at  the  positions  of  the  vibrational  levels  of 
the  02  system.  This  compound  (resonance)  state  decays  by  autodetachnment 
which  accounts10  for  the  vibrational  cross  section  of  O2.  However,  this  state 
can  be  stabilized  in  collisions  with  O2  resulting  in  a  stable  02  .  The 
structures  observed  by  Spence  and  Schulz10  establish  the  validity  of  the 
assumption^1  that  the  three-body  attachment  is  a  two  step  process,  i.e.. 


E 


'2 


E 


(O2)  +  02 


o2  o2 


(8) 


The  structure  in  the  rate  coefficient  has  also  been  observed  by  McCorkle, 
et  al*®  especially  for  electrons  with  temperature  below  0.1  eV  as  shown  in 
Figure  2. 

2.2  THREE-BODY  ATTACHMENT  (M=N2) 

The  attachment  to  02  with  N2  as  the  third  body  has  been  measured  by  many 

workers**’*3’*^’*®  where  the  room  temperature  values  of  a 3  (N2)  are:  5.6  x 

10‘32  cm6/sec  [Ref.  11],  1.5  x  10‘31  [Ref.  16],  and  8.5  x  10"32  [Refs.  13  and 

14].  The  last  value  falls  almost  midway  between  the  other  two  measurements 

and  hence  may  be  used  as  an  average  value  at  room  ’  temperature.  The  most 

recent  measurements*3’*®  which  are  performed  at  high  pressure  (<  300  Torr)  are 

higher  than  the  swarm  data  which  were  performed  at  pressures  much  lower  than 

300  Torr.  These  recent  measurements  indicate  that  at  higher  pressures*3  the 

attachment  rate  data  do  not  reflect  the  three-body  process  and  that  one  has  to 

also  invoke  a  four-body  process*3’*®’22.  These  processes  (see  Eqs.  9  and  10) 

—  * 

follow  those  in  Eqs.  7  and  8  where  the  excited  negative  ion  (02)  is  converted 
into  a  coumpound  of  the  form  [02*  -  N2]  in  collisions  with  N2.  This  compound 
then  is  stablized  through  a  collision  with  02  or  N2,  resulting  in  a 
stable  02  : 


(02)  +  N2  [02  -  N2] 


9 


-\  •*. 


(9) 


[02  "  +  O2  ♦  O2  +  N2  +  O2 


(10a) 


N2  ]  +  N2  >  02  Jr  n2  +  ^2 


(10b) 


The  temperature  dependence  of  aj^)  based  on  the  data  of  References  11, 
16  and  22  are  shown  in  Fig.  2.  McCorkle,  et  al*®  indicate  that  if  they  replot 
the  data  of  Chanin,  et  al*1  by  obtaining  the  mean  energy  from  the  know 
electron  energy  distribution  in  N2,  an  excellent  ageeement  between  such  a  plot 
and  McCorkle  et  al*6  data  is  realized. 

In  the  four-body  process  (see  Eqs.  9  and  10),  it  is  found^  that  02  is 
more  effective  than  N2  in  stabilizing  02  .  The  room  temperature  values  for 
the  four  body  attachment  rate  coefficients  where  02  and  N2  are  the  fourth 
bodies  are7: 


a4  (O2)  =  1.0  x 


-49  9 

10  cm  /sec 


(ID 


On  (N2)  =  6.6  x  10-52  cm^/sec 


(12) 


These  processes  increase  the  three  body  attahcment  in  one  atmosphere  of  air  at 
room  temperature  by  ~  20%. 
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2.3.  THREE-BODY  ATTACHMENT  (M=H20) 

The  three-body  attachment  to  02  with  H20  as  the  third  body  has  been 
measured^3’^  for  electrons  at  room  temperature.  The  coefficients  for  a3 
(H20)  are:  (1.4  +  0.2)  x  10"^9  cm^/sec  [Ref.  23]  and  (1.4  +  0.5)  x  lO"^9 
cm®/sec  [Ref.  24]  which  are  in  excellent  agreement.  However,  no  data  exist 
for  this  attachment  for  electron  temperature  above  400°K. 

2.4  THREE-BODY  ATTACHMENT  IN  DRY  AIR 

The  coefficients  obtained  experimentally  for  the  three-body  and  four-body 
attachment  rates  to  02  with  the  third-body  being  02  and  N2  can  be  utilized  to 
obtain  the  attachment  rate  in  air.  However,  direct  meansurements  of  the 
attachment  rate  in  air  also  exists  and  are  reviewed  in  this  section.  The  room 
temperature  rate  coefficient  for  the  three  body  attachment  in  air,  (air), 
is  (1.1  +  0.07)  x  10'31  cm6/sec  [Ref.  17]  which  is  obtained  in  air  like 
mixture  (4N2:02).  The  electron  life  time  measurements^  in  air  afterglow 
plasma  for  air  pressures  between  50  and  600  Torr  yield  the  following 
expression  (which  may  have  an  uncertainty  of  -  30%) 


(13) 


where  p  is  in  Torr.  Accordingly,  the  life  time  of  thermal  electrons  in  one 
atmosphere  is  1.1  x  10"®  sec.  Similar  results  have  been  obtained  by  Price  and 
VanLint^.  Relation  (13)  implies  that  total  attachment  in  air  at  760  Torr  is 
1.24  x  10“3*  cm®/sec.  If  we  accept  the  average  values  of  83  (02)  and  aj  (N2) 
as  2.4  x  10-3®  cm®/sec  and  8.5  x  10-3^  cm®/sec,  respectively,  the  three-body 
attachment  rate  coefficient  in  air  will  be  1.05  x  10"31  cm6/sec.  If  one 
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adds  ~  20%  increase  due  to  four  body  effects  one  obtains  a  rate  coefficient  of 
1.25  x  lO*^1  cm6/sec  which  is  in  good  agreement  with  the  lifetime 
measurement26  and  is  shown  in  Fig.  3. 

The  temperature  dependence  of  the  three-body  attachment  in  air  has  been 
measured  by  Schneider  and  Brau27  in  the  temperature  range  of  0.45  to  1.2  eV. 
These  measurements  are  shown  on  Fig.  3  along  with  the  calculated  values  from 
the  data  of  Fig.  2  where 

a3  (air)  =  0.04  a3  (02)  +  0.16  a3  (N2)  (14) 

In  this  case  we  have  utilized  the  data  for  N2  by  McCorkle,  et  al16  which  is  in 
excellent  agreement  with  the  data  of  Chanin,  et  al11,  replotted  as  a  function 
of  the  average  electron  energy16.  As  for  02*  the  data  of  Chanin,  et  al11  and 
McCorkle,  et  al16  were  treated  separately  and  the  results  are  shown  on  Fig.  3 
as  a  function  of  the  electron  temperature.  It  is  obvious  from  the  figure  that 
good  agreement  exists  for  the  coefficients  for  Tg<  0.1  .  Above  Tg>  0.1  the 
two  sets  differ  by  as  much  as  a  factor  of  3.  However,  the  set  obtained  from 
the  data  of  McCorkel ,  et  al16  merges  in  a  resonable  manner  with  those  high 
temperature  measurements  of  Schneider  and  Brau27.  It  is  possible  that  if  the 
right  energy  scale  is  utilized  for  aj  (O2)  of  Ref.  11  its  data  will  be  shifted 
towards  lower  temperature  resulting  in  better  agreement  with  the  data  of  Ref. 
16. 

2.5  THREE-BODY  ATTACHMENT  IN  WET  AIR 

There  are  no  temperature  depended  measurements  for  the  three-body 
attachment  to  O2  with  H2O  as  the  third  body.  The  data26,2^  for  thermal 


electrons  at  300°K  are  in  good  agreement  with  a  value  of  1.4  x  10"^ 
cm8/sec.  This  value  indicates  that  H20  is  6  times  more  effective  than  02  as  a 
third-body  in  stabilizing  0^  .  If  we  assume  that  the  temperature  dependence 
of  attachment  with  H20  as  the  third  body  is  similar  to  that  with  02  as  the 
third  body,  then  the  three-body  attachment  rate  with  H20  as  the  third  body 
should  follow  the  data  in  Fig.  2  with  an  increase  in  magnitude  by  a  factor  of 
6.  This  has  been  done  and  the  results  are  shown  in  Fig.  4. 

If  we  utilize  the  room  temperature  value  of  the  three  body  attachment  to 
02  with  H20  as  the  third  body  (1.4  x  10-^cm8/sec)  and  1.25  x  10-31  cm6/sec  in 
dry  air  we  obtain  attachment  rate  for  moist  air.  We  consider  three  values  of 
2%,  4%  and  5%  for  H20  in  air  and  calculate  the  attachment  rate  in  moist  air. 
The  results  are  shown  in  Fig.  5  along  with  the  experimental  results  of  VanLint 
and  Price88. 

3.  THE  ELECTRON  TEMPERATURE 

In  Section  2  we  discussed  the  attachment  rates  in  dry  and  wet  air.  It  is 
obvious  that  these  rates  are  dependent  on  the  electron  temperature  which  has 
to  be  known  as  a  function  of  E/p  for  a  realistic  modelling  of  nuclear  induced 
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lightning.  The  electron  temperature,  defined14-  in  the  swarm  experiments  as 
the  characteristic  energy  is  often  obtained  as  a  function  of  E/p  by  measuring 
Dt/y  .  Where  Dt  is  the  diffusion  coefficient  transverse  to  the  direction  of 
the  electric  field  and  y  is  the  electron  mobility.  The  electron  temperature, 
Te,  in  dry  air  has  been  measured  by  many  investigators8^-33  for  a  wide  range 
of  E/p.  These  measurements  compiled  by  Gallaher  et  al34  are  shown  in  Figures 
[6  a,b]  where  good  agreement  is  apparent,  whenever  more  than  one  measurement 
exists  for  a  given  range  of  E/p.  However,  no  measurement  of  Tg  exists  for 
E/p  below  0.1  V-cm-l-Torr-l .  On  the  other  hand,  Te  has  been  measured  and 
calcualted35-40  in  both  No  and  0?  for  values  of  E/p  far  below  0.1  V  cm-1 


Torr"*.  For  N£  the  data  of  Jory^  and  Crompton  and  Elford^  are  shown  in 
Figure  [6b]  along  with  the  theoretical  values  of  Hake  and  Phelps4^  for  Tg  in 
O2  • 

It  is  interesting  to  note  that  the  electron  temperature  measurements  of 
Crompton  et  al^  in  air  fall  between  the  Te  values  in  N2  and  O2  except  for  the 
data  at  lowest  end  of  E/p.  Furthermore,  the  Te  values  in  N2  and  02  get  closer 
as  E/p  gets  lower  converging  to  the  thermal  value  as  expected.  Assuming  that 
the  diffusion  coefficient  Dt  is  the  same  in  N2  and  O2  one  can  utilize  the  data 
for  Te  in  N2  and  O2  to  obtain  Te  in  air.  This  we  have  done  and  the  result  is 
shown  in  Fig.  [6b]. 

However,  for  most  modelling  purposes  one  is  also  interested  in  the 
electron  tempeature  in  moist  air.  Unfortunately  no  such  measurement  exists. 
On  the  other  hand,  measurements  and  calculations  of  Te  in  pure  water  vapor 
exists.  This  data  is  shown  in  Fig.  [6a]  for  the  sake  of  comparison  with  Tg  in 
dry  air  and  is  obtained  from  the  measurements  of  Parr  and  Moruzzi4*  and  that 
of  Wilson,  et  al4^  whose  D^/u  measurements  agree  well  with  the  theoretical 
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values  of  Lowke  and  Parker  •  From  this  figure  one  sees  that  the  electron 
temperature  rises  very  rapidly  for  E/p  between  15  and  30  V  cm"*  Torr"*.  For 
E/p  >  20  V  cm"*  Torr"*  the  electron  temperature  in  H2O  is  slightly  higher  than 
that  in  dry  air.  However,  for  ^  <  10  V  cm”*  Torr"*,  the  electron  temperature 
in  H2O  is  much  lower  than  that  in  dry  air.  In  this  region 
(p-  <  10  V  cm"*  Torr"*)  the  electron  temperature  in  H2O  is  thermal 
(i.e.  ~  300°k)  while  the  electron  temperature  in  dry  air  is  much  higher  (see 
Fig.  6a).  Between  E/p  *  1.0  to  10.0  V  cm"*  Torr"*  the  electron  temperature  in 
air  is  12  to  40  times  higher  than  that  in  water  vapor.  Therefore,  in  moist 
air,  the  electron  temperature  will  be  lower  than  in  the  case  of  dry  air,  in 


the  E/p  range  just  alluded  to.  If  one  assumes,  for  the  sake  of  argument  that 
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the  diffusion  coefficient  in  ^0  is  the  same  as  in  dry  air  then  for  air  with 
2.5%  of  H2O  the  average  electron  temperature  will  be  lowered  from  1.3  eV  for 
dry  air  to  0.55  eV  which  implies  that  water  vapor  has  a  large  influence  on  the 
electron  temperature  for  E/p  <  10  V  cm"*  Torr~*. 

4.  ION-ION  RECOMBINATION 

In  air  a  variety  of  positive  and  negative  ions  are  formed  as  a  result  of 
the  ionization  of  air  molecules  and  the  subsequent  chemical  processes.  The 
positive  ions  recombine  with  the  free  electrons  and  the  negative  ions.  The 
reocmbi nations  of  the  positive  and  negative  ions  proceed  through  the  two  and 
three-body  neutralization  processes  i.e., 

A+  +  B"  -*■  A  +  B  (15) 

A+  +  B“  +  M>A  +  B  +  M  (16) 

where  A+  and  B"  denote  atomic,  molecular  and  cluster  ions. 

4.1  TWO  BODY  I ON -I ON  RECOMBINATION  (a2) 

The  two  body  mutual  neutral ization  of  positive  and  negative  ions  has  been 
measured  by  the  SRI  group^  using  the  merging  beam  technique.  Most  of  these 

measurments  were  concerned  with  atomic  and  molecular  systems  of  interest  shown 

in  TABLE  II. 


TABLE  II:  TWO  800Y  MUTUAL  NEUTRALIZATION  COEFFICIENT 


In  Units  of 

10’^cm^/sec  (Ref. 

43  &  46)  and  10-^cm^/sec 

(Ref.  44  &  45) 

Reaction 

Ref.  43 

Ref.  46 

Ref.  44  &  45 

0+  +  O’ 

2.7  +  1.3 

N+  +  O’ 

2.6  ±  0.8 

O2  +  0 

1.0  +  0.4 

O2  +  O2 

4.2  t  1.3 

1.0  +  0.1 

N2  +  O2 

1.6  +  0.5 

N0+  +  O2 

5.8  ±  1.0 

nJ  +  NO2 

1.3  +  0.5 

N0+  +  NOi 

5.1  +  1.5 

1.75  +  0.6 

6.4  +  0.7 

O2  +  NO2 

4.1  +  1.3 

N0+  +  NO3 

8.1  +  2.3 

0.34  ±  0.12 

5.7  +  0.6 

O2  +  NO3 

1.3  +  0.4 

H30+  (h20)3 

+  NO3 

5.5 

h30+  (h20)3 

+  no;  +  HN03 

5.7 

However,  the  mutual  neutralizaiton  rate  coefficient,  at  room  temperature, 
given  in  Table  II  by  the  SRI  group4^  are  obtained  by  extrapolation.  The 
actual  measurements  are  at  energies  of  0.15  eV  and  higher.  These  measuremnts 
are  fit  to  a  theoretical  form  predicted  by  Landau-Zener  theory4-*,  and 
extrapolated  to  thermal  energies  and  a  rate  coefficient  is  obtained  from  the 
product  of  the  relative  velocity  of  the  ions  and  the  cross  section  for 
neutralization.  However,  the  afterglow  measurements  of  Smith  et  al44  and 
Smith  and  Church4^,  for  the  reactions  of  N0+  with  NQj  and  NO3  ,  predict 
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neutralization  values  which  are  an  order  of  magnitude  smaller  than  those  of 
Moseley  et  al43.  Eisner  and  Hirsh4®  have  also  measured  these  two  reactions 
and  obtained  [1.75  +  0.6]  x  10"7  cm3/sec  and  (3.4  +  1.2)  x  10'®  cm3/sec  for 
the  neutral ization  of  N0+  with  NO2  and  NO3  ,  respectively.  These  values  are 
lower  from  those  predicted  by  Moseley  et  al43  by  a  factor  of  3  and  30  for  the 
corresponding  neutralizations,  respectively.  The  discrepancy  may  lie  in  the 
extrapolation  methods  of  Ref.  [43]. 

Smith  and  his  colleagues44,4®  have  measured  the  neutralization 
coefficients  for  more  complex  ions  (hydrates),  in  addition  to  those  discussed 
above,  they  find  that  NO3  ions  and  their  clusters,  regardless  of  the  degree 
of  clustering  behave  similarly,  i.e.,  the  mutual  neutralization  coefficients 
of  these  simple  and  complex  ions  are  in  the  range  of  (5-6)  x  10"8  cm3/sec  (see 
Table  II).  This  value  is  in  good  agreement  with  the  measurements  of  Ulwick47 
in  the  upper  atmosphere  (50  -  75  km)  which  give  a  value  of  6.5  x  10-8  cm3/sec. 

The  temperature  dependence  of  the  ion-ion  neutralization  coefficient  has 
the  form  of  (Tg)-8,8  which  is  predicted  theoritical ly48,49,  and  hence  it  can 
be  utilized  to  obtain  the  coefficients  at  temperatures  other  than  at  300°K. 

Theoretical  calculations  of  the  ion-ion  neutralization  coefficients  have 
been  carried  out  by  Olsen48  based  on  a  semiemperical  absorbing  sphere  model. 
The  predicted  values  in  the  majority  of  the  cases,  presented  in  Table  II,  are 
always  below  the  experimentally  obtained  ones.  It  is  interesting  to  note  that 
in  the  case  of  N0+  +  NO3  calculated  value  is  lower  by  a  factor  of  7.  For 
atomic  and  molecular  ions,  however,  experimental  predictions  are  higher  by  a 
factor  of  2  compared  to  the  calculated  values. 

A  middle  ground,  for  modelling  purposes  is  to  use  a  value  of 

1.5  x  IQ'7  (i88)8,8  for  an  neutralization  which  involve  an  atomic  and 
9 

molecular  ion.  However,  whenever  a  triatomic  or  heavier  molecule  is  involved. 


the  measured  value  of  5  x  10"®  cm'Vsec  should  be  utilized.  This  same  value 
should  also  be  utilized  for  complex  ions  (hydrated  ions). 

4.2  THREE  BODY  ION-ION  RECOMBINATION  (<*3) 

The  three  body  ion-ion  recombination  of  various  air  species  ions  of 
interest  have  been  measured  by  McGowan50  at  room  temperature  and  760  Torr. 
These  measurements  are  performed  in  clean  air  where  the  possible  positive  and 
negative  ions  were,  N^,  N3,  N0+,  O4,  O2,  O4,  NOj  . 

The  recombination  rate  coefficient  varied  from  2.28  x  10"®  cm^/sec  to 
2.18  x  10“®cm^/sec  which  imply  a  third  body  rate  coefficient  of  8.5  x  10"2® 
cm®/sec  and  8.10  x  10“2®  cm®/sec,  respectively.  When  air  contained  water 
vapor,  the  ion-ion  recombination  rate  coefficient  varied  from  6.1  x  10“2® 
cm®/sec  to  5.3  x  10“2®  cm®/sec  for  a  relative  humidty  of  38%  and  32%, 
respectively. 

The  value  of  ion-ion  recombination  coefficient  in  one  atmosphere  of  air, 

2.2  x  IO-®  cm^/sec,  obtained  by  McGowan50  is  in  good  agreement  with  that 
measured  by  Sayers51.  However,  a  recent  determination52  of  a3  in  dry  air 
(1  atm,  300°k)  gives  a  value  of  5.0  x  10-7  cm2/sec  which  is  smaller  by  a 
factor  of  4  compared  to  that  measured  by  McGowan50. 

From  the  theoretical  point  of  view,  the  first  attempt  to  calculate 
a3  was  made  by  Thompson5^.  Thompson's  theory  predicts  results  in  good 
agreement  with  experimental  findings  for  pressures  below  one  atmosphere.  The 
success  of  this  theory,  however,  is  largely  fortuitous5^  and  various  attempts 
have  been  made55,5^  to  refine  it.  At  higher  pressures,  on  the  other  hand, 
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Langevin's  Theory  is  more  appropriate  and  has  a  simple  expression  for  the 
recombination  based  on  the  ionic  mobility.  Thus 


Using  the  mobility  dataJJ  for  air  ions,  one  obtains  a3  =  4.7  x  10  cm  /sec  at 
1  atmosphere  which  is  larger  by  a  factor  of  2  compared  to  the  measurements  of 
McGowan50 

5.  IONIZATION  FREQUENCY 

The  ionization  frequency  in  air,  v.  ,  as  a  function  of  E/p  is  another 
important  parameter  for  the  modelling  of  the  nuclear  induced  lightning.  The 
ionization  frequency,  in  general,  is  obtained  from  the  ionization  coefficient 
(Townsend  coefficient)  a  and  the  drift  velocity  Vd  where 

vi  =  a  Vd  (18) 

The  ionization  coefficient  in  air  has  been  measured  by  many 
investigators00,56-59  and  for  a  wide  range  of  E/p.  These  measurements,  (Ref. 
3,  56  to  59)  are  in  reasonable  agreement  especially  those  which  are  for  clean 
air5®-59.  The  old  measurements56,5^  are  not  reliable  at  low  E/p  since  they 
were  contaminated00  with  mercury  which  has  low  ionization  potential.  A  best 
fit  to  the  data  of  Reference  (58)  yields  the  following  expressions60  with  an 
error  of  <  14% 

~  =  8.34  Exp  (-  273.8/E/p)  120  >  £  >  54  (19) 

j  =  16.0  Exp  (-  359/E/p)  1000  >  E/p  >  120  (20) 

The  drift  velocity  in  air,  on  the  other  hand,  can  be  expressed60  as 


54  <  E/p  <;  120 


(21) 


Vd  =  6.0  x  106  +  2.5  x  105  (E/p) 


V.  =  3.38  x  106 
d 


(E/P)l72 


120  <  E/p  <  1000 


(22) 


These  expressions  are  obtained®®  using  the  measured  drift  velocities  in 
nitrogen  and  oxygen.  The  corresponding  ionization  frequency  in  air  is 


v. 


P 


[6.0  x  106  +  2.5  x  105  (E/p)]  8.34  Exp 


for 


54  <  E/p  <  120 


(23) 


for 


=  54.08  x  106  (E/p)1^  Exp(^|) 
120  <  E/p  <  1000 


(24) 


For  the  region  below  =  50  V  cm"*  Torr"*  ,  the  drift  velocity®*’®2  in  air 
fits  the  following  expression 


Vd  =  8.7  x  105  (E/p )0*76 


30  <  E/p  <  50 


(25) 


to  within  10%  of  the  data.  For  E/p  below  54  V  cm'*  Torr'*,  the  ionization 
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coefficient  can  be  expressed  as 


a  _  n  n  r..»  f~  273. 8V 

p  =  12*°  Exp  ( — e7T} 


37  <  —  <  54 
P 


which  fits  the  measured59  data  to  better  than  10%.  Thus  the  ionization 
frequency  in  air  for  E/p  =  37  -  54  V  cm"*  Torr"*  is 


~  =  1.04  x  107  (E/p)0*76  Exp  (--■?-") 


This  may  be  extended  down  to  =  30  V  cm~^  Torr  *  which  is  the  standard 
reduced  breakdown  field  in  air. 

6.  ASSOCIATION  REACTIONS 

In  dry  air  the  primary  ions  generated  are  N2,  02,  N+  and  0+  which  undergo 
certain  association  reactions  that  convert  them  into  heavier  ions,  e.q., 


Nj  +  N2  +  N2  -  nJ  +  N2 


N+  +  N2  +  N2  -  N3  +  N2 


02  +  02  +  02  -*■  04  +  02 


•Tf 


O2  +  ^2  +  ^2 


O2  •  N2  +  N2 


(31) 


N0+  +  N2  +  N2  -  NO+.N2  +  N2 


(32) 


The  dissociative  recombinations  of  these  air  ion  clusters  with  the  plasma 

electrons  proceed  at  a  rapid  rate  compared  to  those  of  the  light  ions  and 

hence  affect  drastically  the  plasma  decay. 

The  association  rate  coefficient  for  reaction  (28)  has  been  measured68-7* 

at  room  temperature  and  as  a  function  of  E/p  where  the  most  recent 

measurement7*  gives  a  value  of  5.0  x  10-^9  cm6/sec.  However,  the  dependence 

of  this  coefficient  on  the  gas  temperature  has  been  investigated  by  Good,  et 

al68  in  the  temperature  range  of  300°k  to  380°k  and  by  Dheandhanoo,  et  'l71, 

in  the  temperature  range  of  120°k  to  480°k.  A  temperature  dependence  of 

Tg“^*°  can  be  inferred7^  from  the  measurement  of  Good  et  al68,  while  the  most 
71  -2  2 

recent  measurement  1  gives  a  Tg  *  dependence  and  it  should  be  preferred 

since  it  is  obtained  over  a  wider  range  of  Tg.  Hence  the  most  current  rate 
coefficient  for  reaction  (28)  is  1.50  x  10~^  (y— cm^/sec  where  Tg  is  in 
eV. 

The  coefficient  for  reaction  (29)  has  been  measured68,69,7**78  and  the 

room  temperature  value,  in  units  of  10-^9  cm6/sec,  is  2.7  [Ref  70],  5.0  [Ref 

68],  1.8  [Ref  73]  and  2.0  [Ref  71].  The  most  recent  measurement7*  also 

provides  the  temperature  dependence  of  the  reaction  which  follows  (Tg)  . 

Hence  the  most  current  coefficient  for  reaction  (29)  is  1.25  x  lO'^y-)2’0  , 

9 

with  Tg  in  eV. 
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The  coefficient  for  reaction  (30)  has  been  measured  at  room 

temperature7**’75  and  as  a  function  of  the  gas  temperature75,75  for  the  range 

of  80°k  to  300°k  (Ref.  75)  and  300°k  to  340°k  (Ref.  76).  The  room  temperature 

value  is  ~  2.5  x  1CT58  cm5/sec  and  the  temperature  dependence  varies  as 

T  ”2*8  .  Thus  the  rate  coefficient  for  reaction  (30)  is  9xl0~35 (T  )~2*8 

g  g 

where  Tg  is  in  eV. 

As  for  reaction  (31),  its  coefficient  and  temperature  dependence,  in  a 

limited  range,  has  been  obtained  most  recently  by  Dheandhanoo,  et  al7*.  The 

36  3  2 

coefficient  is  7.5  x  10  (Tg)  where  Tg  is  in  eV.  However,  for  most 
modelling  purposes  one  can  ignore  O2.N2  »  for  two  reasons.  One,  its 
collisional  breakup  with  N2  (~  2  x  10”^  cm3/sec)75  is  as  rapid  as  its 
formation  (in  one  atmosphere)  and  second  is  its  rapid  switching  reaction77 
with  02  leading  to  the  formation  of  O4  . 

The  last  reaction  in  this  section  for  the  dry  air  clusters  (Eq.  32)  has 
recently  been  meausred71  and  the  coefficient  is  2.68  x  10~38  (Tg)~4*4  where 
Tg  is  in  eV. 

7.  THERMAL  PROCESSES 

The  streamer  propagation  is  initiated  near  a  sharp  metallic  tip  which 
enhances  the  electric  field.  The  ionization  and  heating  of  air  in  this 
region,  which  occurs  as  a  result  of  ohmic  heating,  requires  a  detailed 
treatment  of  the  air  heating  processes.  The  energy  gained  from  the  field  by 
the  electron  is  expended  in  collisions  with  the  air  molecules.  These 
collisions  are  elastic  and  inelastic  in  nature  and  result  in  direct  or 
indirect  heating  of  the  air  species.  The  inelastic  processes  and  their 
contribution  to  the  heating  of  air  molecules  are 


a)  The  Ionization  Energy 

Part  of  this  energy  returns  into  the  gas  through  the  dissociative 


recombination  of  the  molecular  ions.  If  the  ions  are  N2+  and  02  +  ,  the 
contributions  per  recombination  to  air  heating  are  ~  5.8  eV  and  6.96  eV, 
respectively.  On  the  other  hand  if  these  molecules  are  converted  to 
N4+  and  04  +  which  is  highly  probable  when  the  air  is  not  heated  at  the 
beginning,  then  the  contributrions  per  recombination  to  air  heating  are 
~  15.0  eV  and  11.6  eV,  respectively.  This  arises  from  the  fact  that  the 
dissociative  recombination  of  N4  +  and  04  +  result  in  excited  ^  and  O2 
molecules  which  in  turn  are  quenched  by  other  molecules,  resulting  in 
heating  of  the  air  species, 
b)  Excitation  of  Electronic  States 

The  energy  lost  by  electrons  in  the  excitation  of  the  electronic  states  of 

N2  and  O2,  can  be  considered  as  indirect  heating  of  the  air  species, 

because  of  their  rapid  quenching  by  N2  and  O2.  For  example  the 

N2(C,*  v=o)  state  is  quenched  by  N2  and  02  with  a  rate  coefficient  of 

1.12  x  10-11cm3/sec  and  2.9  x  10"10cm3/sec.  In  one  atmosphere  of  air, 

where  most  of  our  concern  is,  these  coefficients  imply  that  the  (C,  0) 

9  -1 

state  is  quenched  with  a  rate  of  ~  1.6  x  10  sec  which  is  ~  50  times 

faster  than  the  radiative  decay  rate  of  the  state.  Hence,  most  of  the 

energy  lost  by  electron  into  the  electronic  states,  of  N2,  especially  the 

triplets  below  12  eV  will  be  converted  into  the  gas  kinetic  energy  of  the 

3 

air  molecules.  However,  the  A  Z  state  may  result  in  the  dissociation  of 


O2.  For  electronic  states,  higher  than  12  eV,  the  most  probable  product 

is  the  dissociation^  of  the  molecule.  However,  the  electron  energy  loss 

to  O2  would  result  mainly  in  the  dissociation  of  the  molecule  and  the 
1  1 

excitation  of  a  A  and  b  Z  states  which  are  metastable  and  are  quenched  in 
time  scales  much  longer  than  ~  1  m  sec  .  However,  they  are  destroyed  by 
electron  deexitation  and  by  detachment  of  02"  .  The  rate  of  energy  loss 


to  various  states  can  be  found  in  Ref.  (60). 
c).  The  Rotational  and  Vibrational  Excitations 

The  electron  energy  loss  to  the  rotational  mode  of  the  molecules  is  an 
inelastic  process  which  indirectly  heats  the  air  moleucles.  The 
rotational  energy  is  rapidly  converted  into  moleuclar  kinetic  energy. 
This  is  due  to  the  fact  that  only  few  collisions  are  needed  for  the 
rotational  energy  relaxation  into  the  thermal  energy. 

However,  the  electron  energy  loss  into  the  vibrational  mode  generates 
an  energy  reservoir  which  does  not  return  rapidly  into  the  thermal  mode. 
Therefore,  it  must  be  accounted  separately  with  a  detailed  consideration 
for  its  coupling  to  the  electron  temperature  and  the  gas  kinetic 
temperature.  This  requires  a  master  equation®®  for  the  vibrational  level 
with  electron  excitation  and  dexcitations ,  vibrational -vibrational  and 
vibrational-translational  trasitions.  The  rate  coefficients  for  the 
excitations  of  the  first  eight  vibrational  levels  of  N2  is  given  in  Table 
III.  These  coefficients  are  obtained®’®1  using  the  measured  cross 
sections0  °  averaged  with  electron  velocity  over  an  electron  Maxwellian 
velocity  distribution.  The  generation  of  excitation  to  higher  levels,  by 
electron  impact,  can  be  obtained  from  the  following  relation 


v+Av 


AV 


(33) 


The  relaxation  of  the  vibrational  energy  through  heavy  particle  collisions 
i.e.,  the  vibrational -translational  transitions,  mainly  occur  in  collision 
of  the  vibrational  levels  with  atomic  oxygen  which  is  the  most  effective 


agent.  The  rate  coefficient05  for  the  deactivation  of  the  v=o  level  by 
oxygen  atom  can  be  expressed  as 


*10 


=  1.1  x 


10"10  Exp  (-69.9/Tg1/3) 


(34) 


where  Tg  is  in  units  of  eV.  The  corresponding  excitation  rate  coefficient 
can  be  obtained  by  detailed  balance.  For  the  vibrational -vibrational 
exchange  process,  the  rate  coefficient®6  is 

10  10  _11  3 /2 

R0’i  =  <ov>  <P0i>  =  3.8  x  10  11  (T  )  1  (35) 

*  9 


Where  the  exchange  probability  <Poi>  is  for  a  harmonic  osillator  and 
symmetric  resonance  exchange.  For  transition  between  different 
vibrational  levels  with  Av  =  ±1  ,  harmonic  oscillator  relations  can  be 
utilized.  Such  a  treatment  not  only  will  describe  appropriately  the 
energy  flow  between  the  vibrational  and  translational  modes  but  also  will 
give  the  molecular  dissociation  through  the  vibrational  ladder. 

7.1  THERMAL  DISSOCIATION  AND  IONIZATION 

In  addition  to  the  ionization  of  air  species  by  electron  impact,  the  air 
molecules  can  be  dissociated  and  ionized  through  thermal  collisions.  The 
thermal  dissociation  rate®7  of  N2  for  reaction  (36) 


N2  +  N2  N  +  N  +  N2 


(36) 


is  1.90  x  10"9  T  Exp  (-9.76/T  )  where  Tn  is  in  units  of  eV  and  is  valid 
g  9  9 

for  the  temperature  range  of  0.5  eV  to  1.3  eV.  The  accuracy  of  this 
coefficient  is  within  a  factor  of  3.  The  thermal  dissociation  of  O2  with  O2 

OO 

and  M2  as  the  collision  partner  have  the  following  rate  coefficients00 
5.8  x  10"9  Tg”0,83  Exp ( -5 . 12/Tg )  and  2.7  x  10‘10  Tg'1,7  Exp  (-5.12/Tg)  , 
respectively,  where  Tg  is  in  units  of  eV  and  their  accuracy  is  within  a  factor 
of  3. 

8.  FINAL  REMARKS 

A  review  of  major  processes  and  their  rate  coefficients  relevant  to  the 
NIL  modelling  is  presented.  Schemes  for  the  calculation  of  the  electron 
energy  flow  into  various  modes  are  also  discussed  along  with  the  basic  air 
heating  processes.  The  uncertainties  in  the  rates  are  discussed  and 
recommended  values  are  given.  Other  chemical  processes  affecting  the  NIL 
modelling  will  be  discussed  in  a  future  report  along  with  recommendations  for 
measurements  of  specific  reactions. 


k*o 


Table  III  -  Electron  impact  excitation  rate  coefficients  of 
eight  ground  state  vibrational  levels  of  N2  (cm^/sec) 


T 

*1 

x2 

*3 

*4 

X5 

X6 

x7 

X3 

0.1 

5.16E-13 

3.30E-16 

1.45E-16 

3.31E-17 

1.34E-17 

2.48E-18 

1.48E-18 

4.49E-20’ 

0.2 

1.77E-11 

3.15E-12 

1.79E-12 

8.31E-13 

4.69E-13 

2.09E-13 

1.37E-13 

1.87E-14 

0.3 

1.42E-10 

5.97E-11 

3.60E-11 

2.04E-11 

1.31E-11 

7.59E-12 

5.15E-12 

1.01E-12 

0.4 

4.58E-10 

2.42E-10 

1.48E-10 

9.16E-11 

6.27E-11 

4.1 IE— 1 1 

2.85E-11 

6.70E-12 

0.5 

9.17E-10 

5.32E-10 

3.27E-10 

2.12E-10 

1.51E-10 

1.05E-10 

7.44E-11 

1.95E-11 

0.6 

1.43E-09 

8.68E-10 

5.32E-10 

3.55E-10 

2.50E-10 

1.90E-10 

1.35E-10 

3.79E-11 

0.7 

1.91 E -09 

1.20E-09 

7.30E-10 

4.96E-10 

3.71E-10 

2.78E-10 

2.00E-10 

5.90E-11 

0.8 

2.33E-09 

1.49E-09 

9.05E-10 

6.23E-10 

4.73E-10 

3.61E-10 

2.61E-10 

8.01E-11 

0.9 

2.67E-09 

1.7  3E-09 

1.05E-09 

7.29E-10 

5.60E-10 

4.34E-10 

3.16E-10 

9.98E-11 

1.0 

2.95E-09 

1.92E-09 

1.16E-09 

8.13E-10 

6.31E-10 

4.94E-10 

3.61E-10 

1.17E-10 

1.1 

3.16E-09 

2.07E-09 

1.25E-09 

8.78E-10 

6.87E-10 

5.42E-10 

3.99E-10 

1.31E-10 

1.2 

3.31E-09 

2.18E-09 

1.31E-09 

9.26E-10 

7.29E-10 

5.79E-10 

4.27E-10 

1.43E-10 

1.3 

3.41E-09 

2.26E-09 

1.35E-09 

9.60E-10 

7.50E-10 

6.07E-10 

4.49E-10 

1.52E-10 

1.4 

3.48E-09 

2.31E-09 

1.38E-09 

9.82E-10 

7.81E-10 

6.27E-10 

4.65E-10 

1.59E-10 

1.5 

3.51E-09 

2.34E-09 

1.39E-09 

9.94E-10 

7.94E-10 

6.40E-10 

4.76E-10 

1.64E-10 

1.6 

3.52E-09 

2.35E-09 

1.40E-09 

9.99E-10 

8.01E-10 

6.47E-10 

4.83E-10 

1.68E-10 

1.7 

3.51E-09 

2.35E-09 

1.39E-09 

9.98E-10 

8.02E-10 

6.50E-10 

4.86E-10 

1.71E-10 

1.8 

3.49E-09 

2.33E-09 

1.38E-09 

9.92E-10 

8.00E-10 

6.50E-10 

4.87E-10 

1.72E-10 

1.9 

3.45E-09 

2.31E-09 

1.37E-09 

9.82E-10 

7.94E-10 

6.57E-10 

4.85E-10 

1.72E-10 

*  4.49E-20  reads  4.49  x  10'20 


TABLE  ill  (Continued) 

Electron  impact  excitation  rate  coefficients  of 
eight  ground  state  vibrational  levels  of  N?  (cm3/sec) 


2.0 

3.41E-09 

2.29E-09 

1 .35  E- 09 

9.7  0E- 10 

7 .85  E- 10 

6.42E-10 

4.82E-10 

1.72E-10 

2.1 

3.36E-09 

2.25E-09 

1.32E-09 

9.55E-10 

7.76E-10 

6.35E-10 

4.78E-10 

1.72E-1Q 

2.2 

3.30E-09 

2.22E-09 

1.30E-09 

9.40E-10 

7.65E-10 

6.26E-10 

4.72E-10 

1.70E-10 

2.3 

3.24E-09 

2.18E-09 

1.28E-09 

9.23E-10 

7.52E-10 

6.17E-10 

4.65E-10 

1.68E-10 

2.4 

3.18E-09 

2.14E-09 

1.25E-09 

9.05E-10 

7.39E-10 

6.07E-10 

4.58E-10 

1.66E-10 

2.5 

3.11E-09 

2.09E-09 

1.22E-09 

8.83E-10 

7.26E-10 

5.97E-10 

4.51E-10 

1.64E-10 

2.6 

3.05E-09 

2.05E-09 

1.20E-09 

8.69E-10 

7.11 E- 1 0 

5.86E-10 

4.43E-10 

1.62E-10 

2.7 

2.98E-09 

2.01E-09 

1.17E-09 

8.50E-10 

7.09E-10 

5.74E-10 

4.35E-10 

1.59E-10 

2.8 

2.92E-09 

1.96E-09 

1.15E-09 

8.32E-10 

6.83E-10 

5.63E-10 

4.27E-10 

1.57E-10 

2.9 

2.85E-09 

1.92E-09 

1.12E-09 

8.14E-10 

6.68E-10 

5.52E-10 

4.18E-10 

1.37E-10 

3.0 

2.79E-09 

1.88E-09 

1.09E-09 

7.95E-10 

6.54E-10 

5.40E-10 

4.1 0E- 10 

1.51E-10 

3.1 

2.73E-09 

1.84E-09 

1.07E-09 

7.78E-10 

6.40E-10 

5.29E-10 

4.02E-10 

1.48E-10 

3.2 

2.67E-09 

1.80E-09 

1.04E-09 

7.60E-10 

6.25E-10 

5.18E-10 

3.93E-10 

1.45E-10 

3.3 

2.61E-09 

1.76E-09 

1.02E-09 

7.43E-10 

6.12E-10 

5.07E-10 

3.85E-10 

1.43E-10 

3.4 

2.55E-09 

1.72E-09 

9.98E-10 

7.25E-10 

5.99E-10 

4.96E-10 

3.77E-10 

1.40E-10 

ELECTRON  TEMPERATURE  (eV) 


The  attachment  rate  coefficient  with  M  =  O2  and  N2  as  a  function  of 
the  electron  temperature. 


ELECTRON  TEMPERATURE  (eV) 


The  attachment  rate  coefficient  in  dry  air  as  a  function  of  the 
electron  temperature. 


ELECTRON  TEMPERATURE  (eV) 


The  attachment  rate  coefficient  with  M  = 
electron  temperature. 


0 t/m  (eV) 


E/P  (V  cm-1  Torr-1) 


Fig.  6b.  The  electron  temperature  in  N£,  O2  and  air  as  a  function  of  E/p 
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